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Abstract 

We prove Lieb-Robinson bounds for systems defined on infinite dimensional Hilbert 
spaces and described by unbounded Hamiltonians. In particular, we consider harmonic 
and certain anharmonic lattice systems. 



1 Introduction 



An important class of systems in statistical mechanics is described by the (an)harmonic lattice 
Hamiltonians, which have a continuous degree of freedom, thought of as a particle trapped in 
a potential, at each site of a lattice. The particles interact by a linear or non-linear force. For 
example, such models are thought to describe the emergence of macroscopic non-equilibrium 
phenomena, such as heat conduction, from many-body Hamiltonian dynamics |24l [2J, the 
understanding of which is one of the long-standing open problems in mathematical statistical 
mechanics [3j. 

Copyright © 2008 by the authors. This article may be reproduced in its entirety for non-commercial 
purposes. 
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In terms of technical difficulty, lattice oscillator models are intermediate between spin 
systems, where the degrees of freedom, each described by a finite-dimensional Hilbert space, 
are labeled by a discrete set, usually a lattice such as Z", on the one hand, and particles 
in continuous space, which necessarily have an infinite-dimensional state space, on the other 
hand. Even in the classical case lattice oscillator systems are significantly more difficult to 
study than spin systems, and also for them more is known than for particle models in the 
continuum. E.g., the existence of the dynamics in the thermodynamics limit was studied by 
Lanford, Lebowitz, and Lieb in [15j . 

In this paper we focus on an essential locality property of the dynamics of quantum 
harmonic and anharmonic lattice models. Since these are non-relativistic models there is 
no a priori bound on the speed of propagation of signals in these systems. In the case of 
quantum spin systems with finite-range interactions, Lieb and Robinson [16] showed that there 
is nevertheless an upper bound on the speed of propagation in the sense that disturbances 
in the system remain confined in a "light" cone up to small corrections that decay at least 
exponentially fast away from the light cone. This is the so-called Lieb-Robinson bound which 
is an upper bound on the speed of propagation. 

In the past few years several generalizations, improvements, and applications of Lieb- 
Robinson type bounds have appeared. This work can be regarded as one further extension, 
going for the first time beyond the realm of quantum spin systems. Here, by quantum spin 
system we mean any quantum system with a finite dimensional Hilbert space of states. For 
example, a quantum spin system over a finite subset A C Z" is described on the Hilbert space 

H A = (g)H x with?4 = <C n * 

where the dimensions 2 < n x < oo are related to the magnitude of the spin at site x G A. 
The algebra of observables for this quantum spin system is then given by 

A A = <g)B(H x )=B(H A ) 

xeA 

where B(7i x ) is the space of bounded operators on 7i x (that is the space of all n x x n x 
matrices). The Hamiltonian of the quantum spin system is usually written in the form 

h a = *P0 

XcA 

where the interaction $ : 2 A — > Aa is such that = 3>(X) G Ax = ® x exB{H x ) for all 

X C A. The time evolution associated with the Hamiltonian H\ is then the one-parameter 
group of automorphisms {r A }i e K defined by 

r t A (A) = jtH AAe -iW A for all A G Aa _ 

For such systems, under appropriate conditions on the interactions $(X) (short-range 
conditions) it was first proved by Lieb and Robinson in [16], that, given A 6 Ax, B 6 Ay, 

\\[T t A (A),B}\\ < Cp||||S|| e -^d(x,Y)-v\t\) ( L1 ) 
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where d(X,Y) = min xe x : y<=-Y \x — y\ and |x| = Y^!j=i \ x j\- The physical interpretation of this 
inequality is straightforward; if two observables A and B are supported in disjoint regions, 
then even after evolving the observable A, apart from exponentially small contributions, their 
supports remain essentially disjoint up to times t < d(X,Y)/v. In other words, this bound 
asserts that the speed of propagation of perturbations in quantum spin systems is bounded. 

In the original proof of the Lieb- Robinson bounds (see [E]), the constant C and the 
velocity v on the right hand side of (II. 1|) depended in a crucial way on N = max xg An x , the 
maximal dimension of the different spin spaces. More recently, new Lieb- Robinson bounds of 
the form (II. ip were derived with a constant C and a velocity of propagation v independent 
of the dimension of the various spin spaces [H], 119] . This new version of the Lieb- Robinson 
bounds allowed for new applications, such as, among other results, a proof of the Lieb-Schutz- 
Mattis theorem in higher dimension, see \12\ I20| . 

It seems natural to ask whether Lieb-Robinson bounds such as (jl.ip can be extended to 
systems defined on infinite dimensional Hilbert spaces, and described by unbounded Hamilto- 
nians. Although the constant C and the velocity v in (jl.ip are independent of the dimension 
of the spin spaces, they depend on the operator norm of the interactions <&(X); for this reason, 
if one deals with unbounded Hamiltonians, the methods developed in [191 CEH1 E] cannot be 
applied directly. Nevertheless, in the present paper we prove that Lieb-Robinson bounds can 
be established for three different types of models with unbounded Hamiltonians, which we 
now present. For the precise statements see Sections [2j O and|H 

First, in Section [21 we consider systems defined on an infinite dimensional Hilbert space by 
Hamilton operators with possibly unbounded on-site terms but bounded interactions between 
sites. In this case, we show that the analysis of |19| goes through with only minor changes, 
and that Lieb-Robinson bounds can be proven in quite a large generality (see Theorem 12. ip . 
A class of interesting examples of this are lattice oscillators coupled by bounded interactions. 
For a finite subset A C Z", one considers the system defined on the Hilbert space Tt\ = 
(£) xeA L 2 (JH,dq x ) by the Hamiltonian 

H = J2p* + V ^ + J2 MO* ~ Qy) 

xeA x,j/GA, \x— 2/1=1 

where p x = —id/dq x , the real function V is such that — A q + V{q) is a self-adjoint operator, 
and <p S L°°(IR) is real valued. Another commonly studied model that satisfies the conditions 
of this result is the so-called quantum rotor Hamiltonian of the form 

Q2 

H = ~ Yl Q#2 + Yl Jx v cos (^ ~ °y + $ 

x x x,y 

where 9 X is the angle associated with the rotor at site x, and J xy are coupling constants 
assumed to vanish whenever \x — y\ exceeds a finite range R. Quantum rotor Hamiltonians 
are used to study a variety of physical situations such as Josephson junction arrays [1], the 
Bose-Hubbard model [22], and crystals consisting of molecules with rotor degrees of freedom 

EH- 

Second, in Section El we consider harmonic lattice systems for which the Hamiltonian 
describes a system of linearly coupled harmonic oscillators situated at the points of a finite 
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subset A cZ", The standard Hamiltonian is of the form 

V 

Rh = Y& + ^ 2 <& + Y Y a j («* - ^) 2 ' 

x [x-i/|=lj'=l 

defined on a finite hypercube in X u , with periodic boundary conditions. In this case, not only 
the on-site terms but also the interactions between sites are given by unbounded operators, 
and the analysis of [19] cannot be applied. As is well-known, the time evolution for harmonic 
systems can be computed explicitly (see Lemma I3.4p , and the derivation of Lieb- Robinson 
bounds (in the form given in Theorem 13. ip reduces to the study of the asymptotic properties 
of certain Fourier sums (see Lemma |3.5() . 

Finally, in Section [H we consider local anharmonic perturbations of the harmonic lattice 
system of the form 

V 

ff = ^ + w 2 ^+ y Y A i fo* - + Y • 

x |x-j)|=lj=l x 

Assuming that the local perturbation V is sufficiently weak (in an appropriate sense), and 
making use of an interpolation argument between the harmonic and the anharmonic time- 
evolution, we derive Lieb- Robinson bounds in Theorem 14.11 

Next, we discuss the classes of observables for which we obtain the Lieb- Robinson bounds 
in each of the three types of models. In the case of quantum spin systems, i.e., the case 
where the Hilbert spaces associated with a lattice site are all finite-dimensional, one proves 
Lieb-Robinson bounds for a pair of arbitrary observables A and B with finite supports (see 
(|l.ip ). It is not clear in general that such a result should be expected when the Hilbert 
spaces are infinite-dimensional and the Hamiltonians unbounded. If the unboundedness in 
the Hamiltonian is restricted to on-site terms while interactions between sites are bounded and 
of sufficiently short range, the standard Lieb-Robinson bound can be derived for arbitrary 
bounded observables. This is explained in Section [2j The novelty of this paper concerns 
harmonic and anharmonic lattice systems which have unbounded interactions of the form 
{lx — Qy) 2 - I n Section [3] and Section S] we prove Lieb-Robinson bounds for Weyl operators. 
The main advantage of working in the Weyl algebra is a consequence of the fact that the class 
of Weyl operators is invariant under the dynamics of the harmonic lattice model, a property 
that is also used in our treatment of anharmonic models. The Lieb-Robinson bounds that we 
obtain for the Weyl operators are sufficient to derive bounds for more general observables, 
such as q x and p x as well as compactly supported smooth bounded functions of q x and p x . 
This is discussed in Section 

Note that locality bounds for harmonic and anharmonic lattice systems have already been 
obtained in the classical setting; while harmonic systems are well-understood, anharmonic 
lattice systems are much more complicated, and a full understanding, even in the classical case, 
has not been reached, yet. In [17], Marchioro, Pellegrinotti, Pulvirenti, and Triolo considered 
anharmonic systems in thermal equilibrium and proved that, after time t, the influence of local 
perturbations becomes negligible at distances larger than i 4//3 . These bounds were recently 
improved in [8] by Butta, Caglioti, Di Ruzza, and Marchioro, who proved that after time t 
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local perturbations of thermal equilibrium are exponentially small in log 2 1 at distances larger 
than t log a t. 

In the quantum mechanical setting, on the other hand, we are only aware of the recent 
work of Buerschaper, who derived, in [7J, Lieb- Robinson type bounds for harmonic lattice 
systems. 



2 Lieb-Robinson Estimates for Hamiltonians with Bounded 
Non-Local Terms 



In this section, we will state and prove our first example of Lieb-Robinson estimates for systems 
with unbounded Hamiltonians. We consider here the dynamics generated by unbounded 
Hamiltonians, assuming, however, the unbounded interactions to be completely local. It 
turns out that, for such systems, locality bounds can be proven in the same generality as 
for quantum spin systems (see Theorem 12.11 below) . Moreover, the proof of this result only 
requires minor modifications with respect to the arguments presented in [19J. 

We first introduce the underlying structure on which our models will be defined. Let T 
be an arbitrary set of sites equipped with a metric d. For T with infinite cardinality, we will 
need to assume that there exists a non-increasing function F : [0, oo) — > (0, oo) for which: 

i) F is uniformly integrable over T, i.e., 

:= supVF(a!(x,y)) < oo, (2.1) 
* er yer 



and 



ii) F satisfies 



C := sup < oo. (2.2, 



Given such a set T and a function F, it is easy to see that for any a > the function 

F a {x) =e- ax F(x), 
also satisfies i) and ii) above with ||-F a || < \\F\\ and C a < C. 

In typical examples, one has that rcZ" for some integer v > 1, and the metric is just 
given by d(x,y) = \x — y\ = Y^j=i \ x j ~ VjY this case, the function F can be chosen as 
F(\x\) = (1 + \x\y v ~ e for any e > 0. 

To each x € T, we will associate a Hilbert space 7i x . Unlike in the setting of quantum spin 
systems, we will not assume that these Hilbert spaces are finite dimensional. For example, 
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in many relevant systems, one considers H x = L 2 (R, dq x ). With any finite subset AcT, the 
Hilbert space of states over A is given by 

H A = (g)?4, 

and the local algebra of observables over A is then defined to be 

where B(H X ) denotes the algebra of bounded linear operators on TC X . 

If Ai C A2 , then there is a natural way of identifying Aa 1 C Aa 2 > and (also in the case of 
infinite T) we may therefore define the algebra of local observables by the inductive limit 

A r = \J A A , 

Acr 

where the union is over all finite subsets A C T; see [U [5] for a general discussion of these 
topics. 

For the locality results we wish to describe, the notion of support of an observable will 
be important. The support of an observable A € Aa is the minimal set X C A for which 
A = A' ® 1 for some A' € Ax = ® x&x B{H X ). 

The result discussed in this section corresponds to bounded perturbations of local self- 
adjoint Hamiltonians. We fix a collection of local operators H loc = {H x } x ^r where each H x 
is a self-adjoint operator over H x . Again, we stress that these operators H x need not be 
bounded. 

In addition, we will consider a general class of bounded perturbations. These are defined 
in terms of an interaction <3?, which is a map from the set of subsets of T to At with the 
property that for each finite set X C V, $(X) G Ax and = <&(X). To obtain our 

bound, we need to impose a growth restriction on the set of interactions $ we consider. For 
any a > 0, denote by B a {T) the set of interactions for which 

||$|| a := sup - 1 Y, II^WI < °°- ( 2 - 3 ) 



Now, for a fixed sequence of local Hamiltonians H = {H x }, as described above, an 
interaction <I> G B a (T), and a finite subset A C T, we will consider self-adjoint Hamiltonians 
of the form 

H A = #i oc + Hf = J2H x +J2 ( 2 - 4 ) 

x<=A XcA 

acting on Ha (with domain given by (£) xe A D(H X ) where D(H X ) C H x denotes the domain of 
H x ). As these operators are self-adjoint, they generate a dynamics, or time evolution, {r t A }, 
which is the one parameter group of automorphisms defined by 

r t A (A) = jtH AAe -itH A for any A 6 A ^ 
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For Hamiltonians of the form (|2.4p . we have a bound analogous to (jl.ip . see Theorem 12.11 
below. 



Before we present this result, we make an observation. It seems intuitively clear that the 
spread of interactions through a system should depend on the surface area of the support of the 
local observables being evolved; not their volume. One can make this explicit by introducing 
the following notation. Denote the surface of a set X, regarded as a subset of A C T, by 

S A (X) = {Z C A : Z n X + and Z D X c ^ 0} . (2.5) 

Here we will use the notation S(X) = Sr(X), and define the ^-boundary of a set X, written 
dg>X, by 

d$X = {x G X : 3Z G S(X) with x G Z and $(Z) / } . 

We have the following result. 

Theorem 2.1. Fix a /oca/ Hamiltonian H^ oc and an interaction G B a (r) for some a > 0. 
Let X and Y" 6e subsets of T. Then, for any A D X U Y and any pair of local observables 
A G Ax and B G Ay, one has that 

\\[t^{A),B]\\ < nA }) m g a (t)D a {X,Y), (2.6) 

where 

e 2||*IUC a |*| _ i ifd{X,Y)>0 



9a{t)-\ e 2\muo a \t\ otherwise, ' (2 ' ?) 



and D a (X,Y) is given by 



D a (X,Y) =min 



53 (<*(*, v)),£ ^ F a (d(x,y)) 



(2- 



The following corollary provides a bound in terms of d(X,Y) = min x€ x,yeY d(x, y), the 
distance between the supports X, Y. 

Corollary 2.2. Under the same assumptions as in Theorem \2.1[ we have 

Mr At a\ Dill / 2 ll^ll W B W W F W ■ rifl x-i ic vii -a n s 

[7f(A),g] < "" mm[|<9$>X| , \d*Y\] e i''> a "J, (2.9) 

Proof of Theorem \2.1[ For any finite Z C T, we introduce the quantity 

m :_ s U p MM, ,2.0) 



z 



and note that C B (Z;0) < 2\\B\\5 Y (Z), where we defined 5 Y {Z) = 1 if YnZ ^ and oy (Z) =0 
if y n Z = 0. A key observation in our proof will be the fact that the dynamics generated by 

H l r + H x = 53^ + 53*(z) 

zGA zcx 
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remains local. More precisely, if we define 

r\ oc {A) = e i( («)4 e -<^+4) for all A £ A A , (2.11) 

we have that for every A £ Ax, tI° c (A) € Ax for every t £ R. This implies, recalling the 
definition (pHO]) . that 

C B {X;t) = sup n . (2.12) 



AeAx 



\A\ 



Consider the function (setting r t (-) = t 4 a (-)) 



f(t) := r t (r^(A)),B 



for A € Ax, B G .4y, and t€R. It is straightforward to verify that 

f{t)=i [rt(HZ))J(t)]-i Y, [rt(r% c (A)),[r t mZ)),B]\. (2.13) 

zes A (x) zes A (X) 

As is discussed in |19[ Appendix A], the first term in the above differential equation is norm 
preserving, and therefore we have the bound 

||/(t)|| < ll/(0)|| + 2|L4|| Yl ^ \\[Ts(HZ)),B]\\ds. (2.14) 

Recalling definition (|2.10p . the above inequality readily implies that 

r\t\ 

C B (X,t) <C B (X,0)+2 Y \M Z )\\ C B (Z,s)ds, (2.15) 

zes(x) Jo 

where we have used (|2.12p . Iterating this inequality, exactly as is done in [19], see also 
[2T] . yields ([12]) with ([22]) and §2JS§. The inequality ([2J|, stated in the corollary, readily 
follows. □ 



In many situations, AcZ" and the bound (|2.9p can be made slightly more explicit (but 
less optimal) by choosing 

F( X ) = (l + \ X \y»-\ andC = 2^^ (i + | 1 x|) ^ +1 . 
In this case we have 

||[t^(A),5]|| < 2- {u+1) \\A\\\\B\\ mm[\d*X\, \d*Y\] e -MV)-2l*II.C|t|), ( 2 .16) 
for all a > 0, with 

||$||a= sup e alx - yl (l + \x-y\Y +1 ||^(^) || < oo . 

x ^ eA X3x,y 

Eq. (|2.16p gives the upper bound 2||<I>|| a C/a for the speed of propagation in these systems. 

One application of the general framework used in Theorem 12.11 concerns systems comprised 
of finite clusters with possibly unbounded interactions within each cluster but only bounded 
interactions between clusters. For such systems, by adjusting T and d(x,y), Theorem 12 . 1 1 still 
applies. 
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3 Harmonic Lattice Systems 



In this section, we present our second example of Lieb-Robinson bounds for systems with 
unbounded Hamiltonians. Let L and v be positive integers. We will consider harmonic 
Hamiltonians defined on cubic subsets Ax = (-L, L] u n 7j u . Specifically, for j = 1, . . . , v and 
real parameters Xj > and u > 0, we will analyze the Hamiltonian 

V 

H h L = H h L {{Xj},u>) = ^p 2 x +u; 2 q 2 x + ^ A, (q x - q x+ej )\ (3.1) 

x€A L i=l 

with periodic boundary conditions (in the sense that q x + ej '■= Qx-(2L-i)e j if z G Ax but 
x + ej Ax), acting in the Hilbert space 

Ka l = (g) L 2 (R,dq x ). (3.2) 

xeA L 

Here {ej}J =1 are the canonical basis vectors in Z u , and since, in most calculations, the values 
of Xj and uj will be fixed, we will simply write for notational convenience. The quantities 
p x and q x , which appear in (13. ip above, are the single site momentum and position operators 
regarded as operators on the full Hilbert space Ha l by setting (we use here units with h = 1) 

p x = t <gi . . . 1 (gi -i— (gi 1 . . . ^ 1 and q x = l®---<g>t®qig)l---®l, (3.3) 
dq 

i.e., these operators act non-trivially only in the x-th factor of Ti.\ L . These operators satisfy 
the canonical commutation relations 

[Px,P y ] = [qx,q y ] = and [q x ,p y ] = i5 X) y, (3.4) 

valid for all x, y G Ax- The Hamiltonian describes a system of coupled harmonic oscillators 
(with mass m = 1/2) sitting at all x G Ax- 

Let A\ L be the algebra of all bounded observables on 7i.A L ■ The time-evolution generated 
by the Hamiltonian (|3.ip is the one-parameter group of automorphisms {t^ ,Al }teR of Aa l , 
defined by 

T?> AL (A)=e UH Ue- UH Z. (3.5) 

As we will regard the length scale L to be fixed, we will suppress the dependence of the 
dynamics on Ax in our notation, by setting t^(.) = r^' L . 



An important class of observables in A\ L are the Weyl operators. For a bounded, complex- 
valued function / : A — > C, we define the Weyl operator W(f) by 

W(f) = e^sA^ 116 /*+P* Im /*) (3.6) 

Clearly, W(f ) is a unitary operator in Aa l such that 

W - 1 (f) = W*(f) = W(-f). 
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Moreover, using the well-known Baker-Campbell-Hausdorff formula 

e A+B = e A e B e -l[A,B] if ^ B]] = [S, [A, B]] = 0, (3.7) 
and the commutation relations (|3.4|) . it follows that Weyl operators satisfy the Weyl relations 
W(f)W(g) = W{g)W(f)e- ilml{f ' 9)] = W(f + g) e -^ lm[{f ' g)] (3.8) 

for any bounded /, g : A — ► C, and that they generate shifts of the position and the momentum 
operator, in the sense that 

W*(f)q x W(f) = q x -lmf x and W*(f) p x W{f) = p x + Re f x . (3.9) 



The main result of this section is a Lieb- Robinson bound for the harmonic time-evolution 
of Weyl operators. 

Theorem 3.1. For any finite X, Y C 1? ' , for all L > such that X, Y C Al, and for any 
functions f and g with supp(/) C X and supp(g) C Y , the estimate 



T?(W(f)),W(g) 



holds for all \i > 0. Here 



< CH/llooNloc £ e -f(«*(».»)-c«.Amax(a, e ^) + i)|t|) (g ^ 



d(x, y) = > min |xj — ?/j + 2Lrjj\ . 



is the distance on the torus. Moreover 



C=[2 + c^l 2 + cZ\ 



(3.11) 



(3.12) 



with c^ x = {lo 2 + 4E" = iA i ) 1/2 '•■ 

Corollary 3.2. Under the same conditions as in Theorem \3.1[ for any < a < 1, one has 

< CH/IUIblU min(|X|,|y|) e ^H^)- c -^ max (f' e(M/2)+1 )l*l) 

(3.13) 



rnW(f)),W(g) 



where 



and 



d(X,Y) = min d(x,y) 

xdX,y£_Y 



-a)\z\ 
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Remark 3.3. i) As we will discuss in Remark \3.6\ (see also Lemma 3.1), both Theorem \3.1\ 
and Corollarv \3.2\ remain valid in the case u = 0. 



ii) If we make the further assumption that the sets X and Y have a minimal separation 
distance, then a stronger, "small time" version of 113. 10\) holds. Specifically, let \i > be 
given, and assume that X andY have been chosen with d(X,Y) > 1 + ^ e • Then for 

any functions f and g with support in X and Y, respectively, one has that 

\[rHW(f)),W(g) 

x£X,y<=Y 

(3-14) 

This bound follows from factoring the t 2 ^ out of \3.43\) , and then completing the argument 
as before. 



< 



Hi) In most applications of the Lieb-Robinson bound it is important to obtain an estimate 
on the group velocity, referred to as the Lieb-Robinson velocity fT8[ \T1\ [T9\ El G21 Gil \2T$. Note 
that we can obtain arbitrarily fast exponential decay in space at the cost of a worse estimate 
for the Lieb-Robinson velocity: 

^) = c^ A max(-,e^ 2 ) +1 ). (3.15) 
A* 

The optimal Lieb-Robinson velocity in the above estimates is obtained by choosing \i = fj,Q, the 
solution of 

2 = e 0*/2)+i . 
I* 

Clearly, 1/2 < fia < 1. This gives the following bound for the Lieb-Robinson velocity in the 
harmonic lattice: Vh(fi ) = 2c LUj \/ fiQ < 4c Wi a- 



Theorem 13. II follows from Lemma [3.41 and Lemma 13.51 both proven below. In Lemma 13.41 
we derive an explicit formula for the time evolution of a Weyl operator. This allows us to 
bound the norm on the l.h.s. of (|3.1U|) by certain Fourier sums which we then estimate in 
Lemma 13.51 

For bounded functions /, g : Al — > C, we define the convolution (/ * g) : Al — > C by 

{f*g)x = fy9x-y, (3.16) 

yeA L 

for any x € (if (x—y) ^ Al, then we define g x -y through the periodic boundary conditions). 

Lemma 3.4. Let L be a positive integer and consider a bounded function f : Al — > C. Then 
the harmonic evolution of the Weyl operator W( f) is the Weyl operator given by 



r t h (W(f)) = W(f t ), ft = f*h^> + f*h%. (3.17) 
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Here the even functions hi t and hi t are given by 



and 



where 



and 



— Im 

2 



|A 



+ Re 



J ik-x-2i-y(k)t 



4%) 



-Im 



D ik-x—2i'y(k)t 



l(k;u),{\j}) = 7(fc) 



+ 4^Aj sin 2 (fc j /2). 
i=i 



(3.18) 



(3.19) 



(3.20) 



The proof of Lemma 13.41 is given in Section 13. 1L 



Lemma 3.5. Suppose that the functions h lt ,h 2t : — ► C are defined as in i3.18\) , h3.19\) . 
Then 



\ht]M)\<(l + \^ /2 + \cZ^ 



2 p (m/2) + 



/or m = 1,2, a// // > 0, i £ M, and x € A^. #ere we defined c^a = (w 2 + 4 5Zj=i ^j) 1//2 an< ^ 
|z| = YTj=i \ x i\- Note that the bounds are uniform in L. 

The proof of Lemma 13.51 can be found in Section 13.21 Using these two lemmas, we can 
complete the proof of Theorem 13.11 



Proof of Theorem \3.1[ Let / and g be functions supported in disjoint sets X and Y, respec- 
tively, with separation distance d(X, Y) > 0. Let L > be large enough so that XUY C Al- 
With Lemma 13.41 and the Weyl relations (|3.8p , it is clear that 

~T t h (W(f)),W(g)} = W(f t ) W(g) (l - e' ilia ^ 

Using the above formula, it follows that 

'r t h (W(f)),W(g)] < | Im [<<?,/*)] | < 



{gJ*h[ L } + f*hf>) 



Expanding the first term, we find that 



(3.21) 



(3.22) 



yeA L 
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and therefore the bound 



< 



oo Hslloo ^2 \ h iJ( x ~y) 

xex,yeY 



(3.23) 



1 



1 



oo y oo 



yv e -Ai(d(*,»)-c„,A max(| , e W 2 )+i) |t|) 



follows from Lemma [3.51 A similar analysis applies to the second term on the r.h.s. of (|3.2ip . 
yielding (fBTTDJ) . □ 



3.1 Harmonic Evolution of Weyl Operators 

The goal of this section is to prove Lemma 13,41 To this end, we diagonalize the harmonic 
Hamiltonian H 1 } by introducing Fourier space operators. Consider the set (recall that Al = 
(-L, L] u n W) 

A i = { X : 1 6 Ai } ' 

Then it is clear that C (— it, 7r] v and |A|J = {2L) U = \Al\. For each k G A£, we introduce 
the operators, 

Qk = and P k = -= == e~ ik - x p x . (3.24) 

V\ A L\ xeAL V\ a l\ x&Al 

One may easily calculate that Q£ = Q- k (similarly, Pi = P- k ) for ah k G A* L . Here we 
have adopted the convention that for k = (k±, . . . , k u ) G A£, —k is defined to be the element 
of A^ whose components are given by 



-kj, if \kj \ < it, 
7r, otherwise. 



This is reasonable as e l7rx = e l7rx for all integers x. These operators satisfy the following 
commutation relations 

[Qk,Qk'\ = [Pk,Pk>] = and [Qk,Pk'\ = ih-k', (3-25) 
for any k, k' G A^- Furthermore, for any x G Al, 

q x = -±= e ik x Qk and p x = — != £ e ik x P k . (3.26) 

With the above relations, it is easy to check that the harmonic Hamiltonian (|3.ip can be 
rewritten as 

H i = 2 PfcP " fc + 7 2 (^)QfeQ-fe • (3.27) 

keA* 
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where we introduced the notation 

7 (fc) = -f(k;{\j},uj) 



\ 



uj- 



+ 4^Aj sin 2 (fcj/2). 
i=i 



(3.28) 



Observe that 7(/c) is independent of sign changes in any component of A;. 



Since we have assumed that uj > 0, we have that ^(k) > to > 0, and therefore, we may 
diagonalize the Hamiltonian by setting 



b k 



Pk-i\l^w-Qk and b* k 



1 P- fc + n/^Q- fe 



V^W) V 2 y^fc) 

In fact, as a result of this definition, we find that for k, k' G 

[h,h'] = [bh b t>] = and [h,b* k ,] = 5 kjk >, 
and moreover, for each k G A 



(3.29) 



(3.30) 



L ' 



Qk 



(b k - b*_ k ) and P k 



l(k) 



Inserting the above into (I3.27p . we have that 



(6, + . (3.31) 



(3.32) 



From this representation of the Hamiltonian H^, we obtain immediately the Heisenberg evo- 
lution of the operators b k and b* k . In fact, from the commutation relations (|3.30p . it follows 
that 

r t h (b k ) = e- 2i ^ k)t b k and T j*(b* k ) = e 2i ^ k)t b* k (3.33) 

for all t G R. 

To compute the evolution of the operators p x and g x , for x G A/,, we express them in 
terms of b k and bt. We find 



Sk-x 



P.i 



L= y e «-p fe = - 7 L= y 



1 v - e 



£ 4ttt ( 6 * - 



(3.34) 



e «- ( bk + 6 * ) 



Therefore 



Sk-x 



' y 



^ e -2i 7 (fc)t &fc _ e 2ij(k)t b *_ 



\ " ^ / ik-x-2if(k)t 7 _ -ik-x+2i-y(k)t i *\ 
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and 

V 2 \ A L\ keAl 

From <^29\\ and (|3T26l) . it follows that 

„• ik-x—2i-y(k)t I i 

A c -ik-x+2i^(k)t I i 

E 7 = 77T ^77T E e lk ' V Py + < V^) E ^% 



2|a z i ,^ y^y i vw) 



feeA* V 'V V \V 'V / ygAi, yGA z 

which implies 



ik-(x—y)—2i-~f(k)t 



Analogously, we find 

2/eA £ 1 L| feeA£ yeA L 1 L| keA* L 

It is then easy to check that 

T t E QxRe /x + Px Im /s I = E 9a;Re +^ Im (/t)x 
with 

and where /i^ and /i^ are defined as in (|3.18|) . (|3.19|) . This proves (|3.17p . 

Remark 3.6. If we consider the Hamiltonian h3. 1\) with to = 0, then we can easily obtain 
analogous formulas for the time evolution of Weyl operators. In fact, if u = 0, we can still 
define operators Pk,Qk is in {3.2J$ and, for every k € A2\{0} ; operators bk and b* k exactly 
as in $3.30\) . In terms of these operators, the Hamiltonian (3. 1\) can be expressed, in the case 
uj = 0, as 

tf£(u; = 0)= P 2 + Yl l{k) {2b%b k + I) . 
keki\{0} 

Since Pq commutes with b^, bt } for all k ^ 0, we obtain (using the commutation relation 13. 30\) 
and t3725\) ) that 

rt(h) = e- 2 ^b k , t*(PD = e^bt 
r^Po) = Po, and T t h (Q ) = Q + 2tP . 
From these formulae, we find that, in the case u = 0, 



r t h (W(f)) = W^f* hfl + / * h { % ) , 
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with 



and where 



and 



|Al| 



~hf}{x) = l -hn 



keA* L \{k } 



+ Re 



J \ " ik-x-2i^{k)t 



|Al| 



fcgA*\{fc } 



ra £ (**>-- 

1 ^ fceA*\{fc } v 



ik-x—2i-*f(k)t 



ik-x—2i-*f(k)t 



(3.35) 



(3.36) 



3.2 Estimates on Fourier Sums. Proof of Lemma 13.51 



The goal of this section is to prove Lemma 13.51 For x E Ax,, let 

i k-x — 2i ■y(k) t 



Hf\t,x) = Re 



— Y 



H { ^(t,x) =Im|^-| 7(ife) e ifc -*- 2i 7(*)f 



j fc-rr — 2 i 7(fc) t 



|A,| fe ^* 7(*) 



(3.37) 



Since /i^(x) = H ( °\t,x) + (i/2) (iT^ (i, zr) + H^fax)) and 4? 0*0 = (i/2) (H^ (t, x) - 



Hi (t,x)), Lemma 13.51 follows from the following exponential estimates on (t,x). 



Lemma 3.7. Suppose that iT^ (i, x), form = —1,1,0, is defined as in (3.37), with j(k) = 
{uj 2 + 4J2 u j=1 Xj sin 2 (Ay /2)) 1 / 2 , and uj > 0. Then we have 

\H^\t,x)\ < e -K^I- c -^ max (f ^ ( " /2)+1 )l*0 

\H^\t,x)\ < c^et e -/«(H-^.A»«(f .^ +1 )l*0 (3.38) 



/or all fi> 0, x e Ax, i € M, and L > 0. ifere c w>A = (w 2 + 4 ^Li Aj) 1/2 . 
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Proof of Lemma \3. 7\ We first prove (|3.38p for m = 0. Since m = throughout this proof, 
and also L is fixed, we will use here the shorthand notation H(t,x) for H^\t,x). We start 
by expanding the exponent e _2 * 7 ( fc )*; 



1 



-2it-f{k)) 7 



n>0 v 7 1 ^ fceA* 

I 1 Nfl ,)nj.2n+l i 

f 2im Z. (2n+l)! |A L | ^ 6 T !M ' 



n>0 



The second term vanishes because 7(— fc) = 7(fc). As for the first term we expand the exponent 
7 2n (/c). We find 



H(t,x) = 



(-l) n 4 n t 2n 



E 



ni 



(2n)\ ^ molmi! . . . m^' 

n>0 ra(),m 1 m„>0 

moH \-m v =n 



2mo 



V 1 

x Yl (4Aj — ^ e ifc ^ sin 2 " 1 ; (Ay /2) 



£=-L+l,...L 



Next we note that, for — L < Xj < L, 



1 £ e^aii^'fc/V) =0 



2L 



-L+1,...L 



(3.39) 



(3.40) 



if \xj\ > rrij. This follows from the orthogonality relation 



- E 

2L ^ 



'ikx 



-L+1,...L 



ii x € A/, , and from the observation that 



e «j*i sin 2m ^(A; i /2) = e* 



(1 — cos kj 



— E 



2 m J 

m,j\ (— 1) 



£=0 



E 

p=0 



# — VP. 



J(xj+2p— i)kj 



(3.41) 



Since — <—£<2p — £<£< rrij, we obtain (|3,40p . Since moreover 



- E 

2L ^ 



< 1 



-L+X....L 
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for all Xj and rrij, we obtain, from f|3.39[) . 

An,f2n „! v 

n>|x| ni|),m 1 ,...,m„>0 J=l 

moH \-m u =n 

v- (2^, A t) 2w 
tT, (2n)! 



(3.42) 



where we put c W) a = (w 2 + 4^J =1 A.,-) 1 / 2 . The previous inequality implies that 

|iHt,x;|S2^ (2n) , S (2|x|) , e . (3.43) 
Using Stirling formula, we find, for arbitrary fx > and for |x| > \t\c w ,\e^/ 2 ) +1 , 

i^,x)i< e - M d a; i-% A i t i). 

Since, by definition |il(t, x)| < 1 for all x € 7j u and t G R, we obtain immediately that 

W>x)| < c -m(W-^.a-«(|,^+ 1 )I*0 

for arbitrary > 0. 

The case m = 1 is handled analogously. For the case m = — 1 we note that 

H i ~ 1 \t,x) = -2 H ( °\s,x)ds (3.44) 
./o 

and then use the bound already obtained for the case m = 0. □ 

4 Lieb-Robinson Inequalities for Anharmonic Lattice Systems 

In this section we consider perturbations of the harmonic lattice system described by the 
Hamiltonian defined in (13. ip . Specifically, for a cube = {—L,L] U C Z", we consider 
the anharmonic Hamiltonian 

H L = H£+ £ 

(4.1) 

= 1] + ^ 2 d + £ Yl a j - ^+^) 2 + Yl v ^ ■ 

xeA L x€A L j=l x£\ L 

We denote the dynamics generated by Hl on the algebra Aa l by r t L ; that is 

Tf L (A) = e itHL A e~ itHL for A G A\ L . 
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The main result of this section will provide estimates in terms of the function 

Since the distance function d is a metric, we clearly have 

Y, F^d(x,z))F^d{z,y)) < C u F^{d{x,y)) 

zeA L 



with 



C,, 



2" +1 V - 

fi + ui)^ 1 ' 



(4.2) 



(4.3) 



zGA L 



Theorem 4.1. Suppose that V € C 1 (R) is real valued with V' € L 1 (M) such that 

Ky = I dw \ V'(w)\\w\ < oo . 



(4.4) 



Then, for every fi>l, and e > 0, there exists a constant C , such that for every pair of finite 
sets X, Y C T, v and L > such that X, Y C Al, we have 



[r t L (W(f)),W(g)} < C||/||oo|b[[ooe^+ e >[*l £ F^d(x,y)) 

x<=X,y£Y 



(4.5) 



for all bounded functions f,g with supp/ C X and supp<? C Y. Here 

(M+j) 

s>0 

and 

v(fi + e) = v h (fi + e) + 
with Vh{fx + e) defined in 113.15]) . 



C = (2 + c^xe^ + c~\) sup [(1 + sY +1 e- 

CC u Ky 



Corollary 4.2. Analogously to Corollarv \3.°A the theorem implies a bound of the form 

[T t L (w(f)),w( g y 

for all [/,, e > and where 

a = r V" L 

(i + M)^ 1 

and d(X, y) denotes the distance between the supports X and Y. 



^ l + U 



Proof. We are going to interpolate between the time evolution r/" (generated by the Hamilto- 



nian (|4.ip ) and the harmonic time evolution t^' Al generated by (|3. 1[) : to simplify the notation 



we will drop all the L dependence in Hl and and in the dynamics r/" and t^' Al . We start 
by noting that 

[rt (W(f)) , W(g)} = \t 8 U_ s (W(f))) , W(g) 



s=t 
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This leads us to the study of 

'rs(rt s (w(f))),W(g) 



d 

ds~ 



T, 



E V(q z ),rt s (W(f)) 

z£A L 

i E [rs([V(q z ),W(f^ s )]),W(g)} 
zeA L 



,W(g) 



(4.6) 



where we used Lemma 13.41 to compute the harmonic evolution of the Weyl operator W(f), 
and the shorthand notation 

= +7*4? (4-7) 

to denote the harmonic evolution of the wave function /. Using (|3.9p . we easily obtain that 



W(Qz),W{f t - B )] = W(f t - a ) (W*(f t - a )V(q z )W(f t - a ) - V{q z )) 
= W{f t - S ) (V(q z - Im f t _ a (z)) - V(q z )) . 



Inserting the last equation in (j4.6f) we find 

£ [r s (r t h _ s (w(f))),W(g) 

= ? 



E \ T » (W(fts) (V{q z - Tmf t - a (z)) - V(q z ))) , W(g)] 
■eA L 

E [ T * { T ts (W(f)j) > W(g)} T S (V(q z - Im/ t _ s (z)) - V(q z )) 



(4.8) 



z&A L 



+ i E r * ( r *- S WW)) ( V (4* - ^ft-s(z)) ~ V{q z )) , W(g)} . 
zeA L 



Next, we define a unitary evolution U{s\t) by 



i—U(s;T) = C(s)U(s;t), and W(r; r) = 1 
as 

with the time-dependent generator 

C(s) = E t 8 (V(q z - 1mft- a (z)) - V{q z )) . 

zeA L 

(Here t > is a fixed parameter). Then, by (|4.8p . we have 

^ [T s (rt s (W(/))),W(5)]w( S ;0) 

= i E ^ (W(/t-.)) [r. (V(q z - lmft-.(z)) - V{q z )) , W{g)\ U(s; 0) 



20 



which implies that 

r t (W(f)),W(g)]u(t;0) 

= [r t h (W(f)),W(g)_ 

V f dsT s (W(f t . s )) [T s (V(q z -Imf t „ s (z))-V(q z )),W(g)]U(s;0) 
STa Jo 



zeA L 1 



Next, we expand 

(V(q z -lmf t . s (z))-V(q z ))= - Tmf t - a (z) f dr V'(q z - r 1mf t -.(z)) 

Jo 

= - hnft-s(z) J dr J dwV'{w)e iw ^- rlraft ~ siz)) 



(4.9) 



where the Fourier transform V is defined as 



V'{w) 



dq 



V'(q)e 



iq-w 



From (|4.9p we obtain 



r t (W(f)),W(g) = r t h (W(f)),W(g) W(0;t) 



- i V f dslmfts(z) f dr [ dwV' (w) e~ iwrlmft ~ s{z) 
zeA L Jo Jo J 

XT s (W(f t - s ))[T s [e iw «*),W(g)]U(s;t). 
Taking the norm, using the unitarity of U(s;t)) and assuming t > for convenience, we obtain 

(W(f)) ,W(g)]\\<\\[r t h (W(f)) ,W(g) 



n 



+ T f t ds\lmf t - a (z)\ [dw\V'(w)\ [r s {e lw ^),W(g)] 
z7A L J ° J 



(4.10) 



For any e > 0, it is clear from (|3.23p that we have 



r t h (W(f)),W(g)\ < (2 + c W)A e^+c- 1 A )||/|| 0O |H| 0O e^+ e H^+ e )* £ e -&*+*M*,v) 

x£X,y£Y 

x£X,y£Y 



where we have set v = (fi + e)v/ l ( / u + e). Similarly, the bound 

\lmft-.(z)\ < Cll/llooe^*-) 



(4.11) 
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follows from an argument as in (|3.23[) . for all < s < t. Plugging these observations into 
(ETTOj) . we find that 



Tt(W(f)),W(g) 
< C 



x<=X,yeY 



zGA L x&X 

Iterating this inequality m times we obtain 
Tt(W(f)),W(g)]\\ 



j dw\V'(w)\ 






Jo 



x£X,y£Y 



<c il/n 

+ C||/|| 00 |b|| 00 e 5i 



E E 



(Ct) n / ^ 



111 



\{J dw 3 \ W] \\V'{ W] )\ 



x£X,y£Y n=l 

x Y F^(d(x,z 1 ))F^(d(z 1 ,z 2 ))...F fM (d(z n ,y)) 
zi,...,z n eA L 

+ c m+i n/iioo e (n / dw j\ w j\\v'{wj)\) f d S1 r d S2 . 



xGX \j=l 



JO 



ds m +i 



x Y F f j,(d{x,zx))F fl (d(zi,z 2 \) ...F^(d{z m ,z m +i)) 

2l,-,2m + l£Ai 



(4.12) 



Using (g2D, we find that 



Y F^d(x,z 1 ))F^d(z 1 ,z 2 ))...F^d(z n ,y)) < C n v F^d(x,y)). 

zi,...,z n eA L 

As for the error term in (|4.12p . we can use the a-priori bound || [T Sm+1 (^e lWrn+iqZm + 1 ), VF((7)] || < 2 
to obtain 

2\\f\\ 00 e^\\V'\\ 1 Ct {C ^ C ^ m Y E F M^*m + l)) 



(m + 1)! 



x£X z m+1 GA L 



< 2 



o C * 111/' 



(m + 1)! 



-\x\ YW)- 

Z&1 V 
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From (|4.12p . we now conclude that 



n(W(f)),W(g)\ < CH/llooll^lUe^ E F^d(x,y))J2 



(CK V C u t) n 



xex,yeY 
+ 2\\f\\ 00 e vt \\V'\\ 1 Ct 



n>0 



(m + 1 ! ^— ' 



<c||/|| 0O y| ooe ( e+o ^^) t E F M^y)) 

x<=X,y£Y 

+ 2 H/IU e « l^lkCt ^^ f |X| £ . 

(m + 1 ! ^-^ 



Since this is true for every m > 0, and since the last term converges to zero as m — » oo, the 
theorem follows. □ 



Remark 4.3. Exactly the same proof yields the Lie}) -Robinson bounds ( f^.5| j for the Hamilto- 
nian 

V 

H L =Y J pl + " 2 al+ EE A ; (& - + E • 

x-eA L ^gA l j=i xeA L 

Moreover, one can see from the proof that the on-site nature of the anharmonic perturbation 
does not play an important role here. For example the same technique can be used to establish 
Lieb-Robinson bounds for the dynamics generated by the Hamiltonian 



Hl= E Px + ^ & + EE A J fe-fe+e,) 2 + E E ( Vl ( qx ~ qx + e j"> + V2 ( Px ~ P»-Hy)) 
x£A L x£A L j=l x&K L j=l 



if both V\ and V2 satisfy the assumption \4-A 



5 Discussion 



5.1 Other Observables 

Theorems 13.11 and 14.11 give a Lieb-Robinson bound for Weyl operators of the form 

\\[Tt(W(f)),W(g)]\\ < CH/IUIblUe-^^)-^) (5.1) 

for / and g supported on finite subsets X and Y of the lattice, where tj is the dynamics of a 
harmonic or anharmonic lattice system that satisfies the conditions of these theorems. From 
(|5,ip one can of course immediately obtain a bound for observables A and B that are finite 
linear combinations of Weyl operators by a simple application of the triangle inequality. Two 
other classes of observables for which we can obtain useful bounds are worth mentioning. 

Note that for every / : X — » C, W(f) = e lb ^\ with a self-adjoint operator b(f) acting on 
Ti-x (|3.6p . such that b(sf) = sb(f) for every s€R. Let A, £> € L 1 (IR) be two functions such 
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that sA(s) and sB(s) are also in L 1 (M). Then, it is straightforward to derive a Lieb- Robinson 
bound for the observables A(b(f)) and B(b(g)) defined by 

A(b{f)) = J dsA(s)W(sf), B{b{g)) = f dsB(s)W(sg). (5.2) 

The result is 

\\[r t (AW))),B(b(g))]\\ < C\\f\U\g\Ue~^ x ^~^ f ds\sA{s)\ f ds\sB{s)\ (5.3) 

By taking derivatives, we can also obtain a Lieb-Robinson bound for the unbounded 
observables b(f) and b(g) (e.g., q x and p x ). Because b(f) and b(g) are unbounded we apply 
the Lieb-Robinson bound first on a common dense domain of analytic vectors (see [5l Lemma 
5.2.12]), and find that the commutator [Tt(b(f)),b(g)] has a bounded extension with the 
following norm bound 

\\[rtW))M9)}\\ ^Cll/Hooll^lloce-^W^I). (5.4) 



5.2 Exponential Clustering Theorem 

For a large class of quantum spin systems it was recently proven that a non- vanishing spectral 
gap implies exponential decay of spatial correlations in the ground state |18|, [T4"l 121 j . Such 
a result is often referred to as the Exponential Clustering Theorem. The locality property 
of the dynamics provided by a Lieb-Robinson bound is one of the main ingredients in the 
proof of this result. In the harmonic case, the clustering properties of the exact ground state 
can be explicitly analyzed [HI [23] , and indeed one finds exponential decay whenever there is 
a non-vanishing gap. For the harmonic systems considered here, the gap is non-vanishing iff 
lo > 0. The results of this paper can be used to prove an exponential clustering theorem for 
the class of anharmonic lattice systems we consider here. In fact, following the method of [21] 
(see also |18t I14j). the only additional estimate needed is the following short-time bound. 

Lemma 5.1. Let Hi be the Hamiltonian acting on Al = (-L, L] u C 7Lf defined in fl^. and 
Tt the time- evolution generated by Hi. Let f,g : — > R with supp/ C X, suppg C Y, and 
X n Y = 0. Then there exists a constant C = C(X,u, Ky) < oo such that 

\\[n(W(f)),W(g)]\\ < C \t\ mm(\X\, \Y\) |||/||oo|bl|oo (5.5) 

for all \t\ < to(X,u, Ky)- 

Proof. Let H L m \t,x), for m = 0, ±1, be the Fourier sums defined in (|3.37p . From (|3.43p . we 
obtain that, for arbitrary /i > 0, 

\H(°)(t,x)\ < (2 Cu)jX \t\) {2C "^ hl < |t| eWO+i e -*M-^l4) 
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for all \x\ > 1 and \t\ < e ^/ 2 ^ 1 c x 1 . Since similar estimates hold for and ^ as well, 

I I — I I A,Lu 1 

we find, analogously to f)3.23j) . that, if rf 1 denotes the harmonic time-evolution generated by 
the Hamiltonian (13.11). 



r t h (W(f)),W(g) 



<C|t|||/||oo||0||oo Yl e 

x£X,y£Y 

< C\t\\\f W^WgW^minilXHY] 



fj,(d(x,y)- 



(5.6) 



for all \t\ < e 1 c uJ \ (using the assumption that X PiY = 0), and for a constant C 

depending only on A and u. 

Next we consider the anharmonic time evolution n = t± . From (]4.10p . it follows that 



(W(f)),W(g)]\\<\\[rHW(f)),W(g) 

+ Y fds |Im/ t _ s (z)| J dw\V'{w)\ [j [ 



t, e 



iwq 



),W(g)] 



(5.7) 



Applying (|5.6p to bound the first term, (14.111) and Corollary 14.21 to bound the second term, 
we find 



rt(W(f)),W(g) 



<C\t\ 



\g\\oo min(|X|, \Y\ 



(5.8) 



for a constant C depending only on A, a; and on the constant Ky defined in (|4.4p . and for all 
\t\ sufficiently small (depending on A, u, and Ky). □ 

As a consequence of these considerations one obtains the following theorem. 

Theorem 5.2. Let H be the Hamiltonian of a harmonic or anharmonic lattice model satisfy- 
ing the conditions of Theorem \3.1\ or \4.1\ and suppose H has a unique ground state Q and a 
spectral gap 7 above the ground state. Denote by {■) the expectation in the state O. Then, for 
any functions f and g with supports X and Y in the lattice we have the following estimate: 

\(W(f)W(g)) - (W(f))(W(g))\ < C||/|U \\g\U ||V^|| X min(|X|, \Y\)e~ d ^ Y ^ (5.9) 

where [i > 1 and e > are as in Theorem \4-l\ and £ can be taken to be 

2(/i + e)v(n + e) + 7 



(5.10) 



and where, if we assume d(X, Y) > C is a constant depending only on the lattice. 



It is straightforward to see that the same bound holds for infinite systems if the corre- 
sponding GNS Hamiltonian has a unique ground state and a spectral gap above it, and the 
infinite system is the thermodynamic limit of finite systems that satisfy the conditions of 
Theorem O or gU 
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